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Abstract The intensity of an electromagnetic wave interacting self-consistently with a beam of charged 
particles as in a free electron laser, displays large oscillations due to an aggregate of particles, called the 
macro-particle. In this article, we propose a strategy to stabilize the intensity by re-shaping the macro- 
particle. This strategy involves the study of the linear stability (using the residue method) of selected 
periodic orbits of a mean-field model. As parameters of an additional perturbation are varied, bifurcations 
occur in the system which have drastic effect on the modification of the self-consistent dynamics, and in 
particular, of the macro-particle. We show how to obtain an appropriate tuning of the parameters which 
is able to strongly decrease the oscillations of the intensity without reducing its mean- value. 

PACS. 94.20.wj Wave/particle interactions - 05.45. Gg Control of chaos - 11.10. Ef Hamiltonian approach 



1 Introduction 

The self-consistent interaction between an electromagnetic 
\ wave (or a set of electromagnetic waves) and a beam of 
charged particles is ubiquitous in many branches of physics, 
e.g. accelerator and plasma physics. For instance, it plays 
a crucial role in the Free Electron Laser, which is used 
to generate a tunable, coherent, high power radiations. 
Such devices differ from conventional lasers in using a rel- 
ativistic electron beam as its lasing medium. The physical 
' mechanism responsible for the light emission and ampli- 
fication is the interaction between the beam and a wave, 
which occurs in presence of a magnetostatic periodic field 
generated in an undulator. Due to the effect of the mag- 
netic field, the electrons are forced to follow sinusoidal tra- 
jectories, thus emitting synchrotron radiation. This initial 
seed, termed spontaneous emission, acts as a trap for the 
electrons which in turn amplify it by emitting coherently, 
until the laser effect is reached. 

The coupled evolution of radiation field and N par- 
ticles can be modeled within the framework of a simpli- 
fied Hamiltonian picture [J]. The N + 1 degree of free- 
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dom Hamiltonian displays a kinetic contribution, associ- 
ated with the particles, and a potential term accounting 
for the self-consistent coupling between the particles and 
the field. Hence, direct inter-particle interactions are ne- 
glected, even though an effective coupling is indirectly pro- 
vided because of the interaction with the wave. 

The linear theory predicts QJ, for the amplitude of the 
radiation field, a linear exponential instability and a late 
oscillating saturation. Inspection of the asymptotic phase- 
space suggests that a bunch of particles gets trapped in 
the resonance and forms a clump that evolves as a single 
macro-particle localized in phase space. The untrapped 
particles are almost uniformly distributed between two 
oscillating boundaries, and populate the so-called chaotic 
sea. 

Furthermore, the macro-particle rotates around a well 
defined center and this peculiar dynamics is shown to be 
responsible for the macroscopic oscillations observed for 
the intensity [SKI- It can be therefore hypothesized that 
a significant reduction in the intensity fluctuations can 
be gained by implementing a dedicated control strategy, 
aimed at confining the macro-particle in space. As a side 
remark, note that the size of the macro-particle is di- 
rectly related to the bunching parameter, a quantity of 
paramount importance in FEL context 0. 

For example, a static electric field |4"ll5llfi] can be used 
to increase the average wave power. While the chaotic 
particles are simply accelerated by the external field, the 
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trapped ones are responsible for the amplification of the 
radiation field. 

The dynamics can also be investigated from a topologi- 
cal point of view, by looking at the phase space structures. 
In the framework of a simplified mean field description, 
i.e. the so-called test-particle picture where the particles 
passively interact with a given electromagnetic wave, the 
trajectories of trapped particles correspond to invariant 
tori, whereas unbounded particles evolve in a chaotic re- 
gion of phase-space. Thus, the macro-particle corresponds 
to a dense set of invariant tori. Our strategy is to mod- 
ify the macro-particle dynamics by restoring or destroying 
invariant tori in selected regions of phase space. 

A technique of Hamiltonian control can be used 0|H] 
to reconstruct additional invariant tori around the macro- 
particle, in order to enhance the trapping. A specific per- 
turbation is computed, which guarantees the confinement 
on invariant tori of trajectories characterized by a specific 
energy. 

In this paper, we propose a strategy to stabilize the 
intensity of the wave, by modifying the characteristics of 
the macro-particle. A (generic) one or several-parameter 
family of perturbations is introduced, which allows us to 
modify the topology of phase-space by tuning appropri- 
ately the parameters. The residue method 11 II 121113] is 
used to identify the important local bifurcations happen- 
ing in the system when the parameters are varied, by an 
analysis of linear stability of selected periodic orbits. This 
technique enables to monitor the size, gyration and inter- 
nal structure of the macro-particle. An appropriate tuning 
of the parameters is able to strongly decrease the oscilla- 
tions of the intensity without reducing its mean- value. 

The paper is organized as follows : In Sec. the test- 
particle model is presented. The latter provides a simple 
topological representation of the self-consistent interac- 
tion. In Sec. 03 the residue method is discussed: in par- 
ticular we show how local bifurcations can generate an 
enlargement of the macro-particle. In Sec.Q] we apply the 
above method to re-shape both the size and the internal 
structure of the macro-particle. Finally, in sect El we draw 
our conclusions. 



2 Dynamics of a single particle 

The dynamics of the wave-particle interaction has been de- 
scribed in Ref. [T] by the TV-body Hamiltonian, accounting 
for a kinetic contribution and an interaction term between 
the particles and the radiation field : 

N 2 j—f N 

H N ({0 jlPj }, <j>, I) = J2 f - 2 Jj^ Yl cos + B i)> (!) 

where (8j,Pj) are the conjugate phase and momentum of 
the z-th particle, and (0, I) stand respectively for the con- 
jugate phase and intensity of the radiation field. Since 
4> is a phase, ((f), I) belongs to T x M+ where T is the 
one-dimensional torus, here (6j,pj) belongs to T x R. The 
phase space of the system is then T N+1 x R w x R + . We 
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Figure 1. Normalized intensity I/N from the dynamics of 
Hamiltonian (J, with N = 10000 particles and Hn = 0, Pn = 

io- 7 . 
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Figure 2. Snapshot of the N particles at t = 1000, with N = 
10000 (same parameters as Fig.IQ. The grey points correspond 
to the chaotic particles, the dark ones to the particles in the 
macro-particle. 

notice that there are two conserved quantities : Hn and 
the total momentum Pn = I + J2j Pj ■ We consider the 
dynamics given by Hamiltonian Q on a 2A^-dimensional 
manifold (defined by Hn = and Pn = s where e is 
infinitesimally small). 

Starting from a negligible level (I — e small and pj — 
0), the intensity grows exponentially and eventually at- 
tains a saturated state characterized by large oscillations, 
as depicted in Fig. Concerning the particle dynamics, 
around half of them are trapped by the wave and form the 
so-called macro-particle (see Fig. El- A particle is trapped 
if it travels in the potential well of the wave, after the 
system has reached saturation, which reads : 

3k e Zsuch thatW > t sat , (2fc-l)7r < 9{t)+<j>{t) < (2fc+l)7r, 

(2) 

The remaining particles experience a chaotic motion within 
an oscillating waterbag, termed chaotic sea, which is un- 
bounded in 9 contrary to the macro-particle. 

In order to get a deeper insight into the dynamics, we 
consider the motion of a single particle. For large N, we 
assume that its influence on the wave is negligible, thus 
it can be schematized as a passive particle in an oscillat- 
ing field. The motion of this test-particle is described by 
the following Hamiltonian with one and a half degrees of 
freedom : 
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Figure 3. Fourier transform of the interaction term h(t), as 
obtained from simulations of Hamiltonian Q, after saturation 
has been reached (same parameters in Fig.Ul. 




Figure 4. Poincare section of a test-particle, described by 
Hamiltonian The periodic orbits with r.n. 7 are marked 
by plus (elliptic orbit) and crosses (hyperbolic orbit). 



Hi p {e,p,t) = P --2^j I -^ C o S {d + <t>{t)) (3) 

= P 1-Re{h{t)e ie ), (4) 

where the interaction term h(t) is derived from some sim- 
ulations of the original iV-body Hamiltonian Q. In the 
saturated regime, hit") is mainly periodic (see Fig. 0). A 
refined Fourier analysis shows that it can be written 
as : 



h{t) = -jJ-e 1 ^ « [F + ae luJlt + f3e- iuit ]e im , (5) 

where fl = —0.685 stands for the wave velocity and u>i — 
1.291 for the frequency of the oscillations of the intensity. 
As for the amplitudes, the Fourier analysis provides the 
following values : F = 1.5382 - 0.0156i, a = 0.2696 - 
0.0734i and (3 = 0.1206 + 0.0306i. 

Hamiltonian Q results from a periodic perturbation 
of a pendulum described by the integrable Hamiltonian 
Ho 

v 2 

i/ = y - \F\ cos(6» + fit + fa), 

where F = |F|e^ F . The linear frequency of this pendulum 
is \J\F\ ps 1.240 which is very close to the frequency of 
the forcing w\ . Therefore a chaotic behaviour is expected 
when the perturbation is added even with small values of 
the parameters a and (3. 

The Poincare sections (stroboscopic plot performed at 
frequency Lof) of the test-particle (see Fig.QJ reveals that 
the macro-particle corresponds to a set of invariant tori. 
Conversely, the chaotic sea is filled with seemingly erratic 
trajectories of particles. The rotation of the macro-particle 
and the oscillations of the waterbag can be visualized by 
translating continuously in time the stroboscopic plot of 
the phase space. 

The macro-particle is organized around a central (el- 
liptic) periodic orbit with rotation number (r.n.) 1. The 
period of oscillations of the intensity is the same as the 



one of the macro-particle which indicates the role played 
by this coherent structure in the stabilization of the wave. 

Thus, in the test-particle model, the macro-particle is 
formed by particles which are trapped on two-dimensional 
invariant tori. This picture can be extended to the self- 
consistent model, if one considers the projection of a tra- 
jectory (cj)(t),I(t),{9j(t),pj(t)}j) in the (0,p) plane, each 
time it crosses the hyperplane J^j sm (4> + = 0, i.e. 
dl/dt = 0. From the full trajectory, we follow a given 
particle (an index j) and plot (9j,Pj) each time the full 
trajectory crosses the Poincare section. 

The trapped particles (following definition (J2J) appear 
to be confined to domains of phase-space much smaller 
than the one of the macro-particle (see Fig. . These do- 
mains are similar to the invariant tori of the test-particle 
model, although thicker. It is worth noticing that not only 
these figures have a similar overall layout but there is 
a deeper correspondence in the structure of the macro- 
particle. For instance, both figures show a periodic orbit 
with period 7 at the boundary of the regular region. Since 
we saw that the macro-particle directly influences the os- 
cillations of the wave, the test particle Hamiltonian Q 
serves as a cornerstone of our control strategy which con- 
sists in modifying the structure of the macro-particle in 
order to stabilize the intensity of the wave. This strat- 
egy focuses on restoring or breaking-up invariant tori to 
reshape the macro-particle. In order to act on invariant 
tori, we use the periodic orbits which, as we have seen in 
Figs. 4 and 5 structure the macro-particle. 



3 Residue method 

The topology of phase space is investigated by analysing 
the linear stability of periodic orbits. Information on the 
nature of these orbits (elliptic, hyperbolic or parabolic) 
is provided using, e.g., an indicator like Greene's residue 
[flHTo] . a quantity that enables to monitor local changes of 
stability in a system subjected to an external perturbation 
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Figure 5. Poincare section of Hamiltonian Q, when the parti- 
cles intersect the plane dl(t)/dt = 0. The different trajectories 
are represented by different grey levels. 

Let us consider an autonomous Hamiltonian flow with 
two degrees of freedom which depends on a set of param- 
eters A G R m : 

z = IWH(z;X), 

where z = (p,E,6,t) £ R 4 and JJ = and ^ 

being the two-dimensional identity matrix. In order to 
analyse the linear stability properties of the associated 
periodic orbits, we also consider the tangent flow written 

as 

j t J l (z) =JV 2 H(z;\)J t , 

where J° = I4 and V 2 H is the Hessian matrix (composed 
of second derivatives of H with respect to its canonical 
variables). For a given periodic orbit with period T, the 
linear stability properties are given by the spectrum of the 
monodromy matrix J T . These properties can be synthet- 
ically captured in the definition of Greene's residue : 



4 - tr J T 



R = 



In particular, if R e]0, 1[, the periodic orbit is called el- 
liptic (and is in general stable); if R < or R > 1 it is 
hyperbolic; and if R = and R = 1, it is parabolic while 
higher order expansions give the stability of such periodic 
orbits. 

Since the periodic orbit and its stability depend on 
the set of parameters A, the features of the dynamics 
will change under apposite variations of such parameters. 
Generically, periodic orbits and their (linear or nonlinear) 
stability properties are robust to small changes of parame- 
ters, except at specific values when bifurcations occur. The 
residue method [TTlll2lll3| detects the rare events where 
the linear stability of a given periodic orbit changes thus 
allowing to calculate the appropriate values of the param- 
eters leading to the prescribed behaviour of the dynamics. 
As a consequence, this method can yield reduction as well 
as enhancement of chaos. 



First, we illustrate the method by introducing an ad- 
ditional interaction term which depends on a parameter. 
The latter has to be properly tuned in order to control the 
dynamics. Here, the control term is chosen as 



//',„( 0. /'. / : Ai = IhjH. ,, 1 s - 2Ay ^ cos (26 + 4>(t)), 

(6) 

and it is therefore similar to the original interaction term 
between the charged particle and the wave. 

Alternatively, other types of perturbations could be 
selected, our choice being solely motivated by didactic 
reasons. For A = (which corresponds to the original 
Hamiltonian Hi p ), we consider two Birkhoff periodic or- 
bits which originates from the breakup of invariant tori 
with rational ratio of the integrable case (given by the 
nonlinear pendulum described by Hq in the previous sec- 
tion). These two periodic orbits have the same action but 
different angles in the integrable case Hq and have the 
same r.n. on the Poincare section; one is elliptic O e and 
one is hyperbolic Oh (see Fig.^Jl. Let us recall that the ro- 
tation number (or winding number) of a periodic orbit is 
the number of times it crosses the Poincare section before 
closing back onto itself. 

We call R e and Rh the residues of these orbits. We 
have i? e (0) > and Rh(0) < 0. We then modify contin- 
uously the parameter A starting from 0. For each value 
of A, we follow continously the periodic orbit under con- 
sideration as well as its linear stability property indicated 
by its residue. We plot the values of the residues as a 
function of the parameter A in Fig. We notice that at 
A = A c m —0.0370 we have : 



Re{Xc) = Rh{\c) = 0. 



(7) 



The bifurcation is associated with the creation of an 
invariant torus [E]. This diagnostic is confirmed by the 
Poincare section (see Fig. EJ of the controlled Hamilto- 
nian ©, at A = A c : The elliptic islands with r.n. 7 have 
been replaced by a set of invariant tori in its neighborhood, 
leading to an enlargement of the macro-particle. Note that 
elliptic islands with r.n. 5 are now present around the reg- 
ular core (these orbits were present in the case A = 0, 
but were both hyperbolic in the chaotic sea). The asso- 
ciated couple of elliptic/hyperbolic orbits can be treated 
similarly as those of r.n. 7 in order to gain further enlarge- 
ment of the macro-particle : we modify the value of the 
parameter A of Hamiltonian Q around A = A c , until the 
condition Q holds for the residues of the considered or- 
bits. Such a condition is verified when A = A' c ~ —0.0746, 
and leads to the apparition of a new set of invariant tori 
(not shown here); again it appears that elliptic islands 
with r.n. 4 surround the macro-particle. The process of in- 
creasing the size of the macro-particle can be iterated one 
step further. We modify the parameter A around A = X' c 
such that condition {JJ holds for the periodic orbits with 
r.n. 4. The residue method predicts the formation of an in- 
variant torus at A" w —0.1321. Inspection of phase-space 
(see Fig.lHJ corroborates our findings, resulting in an ad- 
ditional extension in size of the macro-particle. 
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Figure 6. Residue of the elliptic (bold curve) and hyperbolic 
points with r.n. 7 of Hamiltonian I0. 
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Figure 7. Poincare section of Hamiltonian JUJ with A = A c ~ 
—0.0370. The periodic orbits with r.n. 5 are marked by plus 
(elliptic orbit) and crosses (hyperbolic orbit). 
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Figure 8. Poincare section of Hamiltonian with A = A" ~ 
-0.1321. 



The test-particle approach is validated using the full 
self-consistent Hamiltonian involving N particles (with 
N ^> 1) interacting with the wave. The control is nat- 
urally introduced in the self-consistent dynamics as : 



H C N {{0j,Pj}, 0, 1] A) = HN({8j,Pj}, <f>, I) 



A 





-0.037 


-0.0746 


-0.1321 


R 


0.7 


0.7 


0.8 


0.9 


Rgyr 


0.17 


0.15 


0.11 


0.03 


(I) 


0.62 


0.66 


0.79 


0.83 


Al 


0.73 


0.66 


0.66 


0.44 



Table 1. Characteristics of the macro-particle and the wave 
for Hamiltonian 10 • 



- 2A V^E C0S ( 2 ^ + ^)' ( 8 ) 

j 

where Hm is given by Eq. . The behaviour of the system 
is investigated for the three critical values where bifurca- 
tions happen in the test-particle model A c , A^,, and A". For 
all these cases, simulations based on Eq. (jHJ display a qual- 
itative behaviour similar to the one obtainde for the orig- 
inal unperturbed dynamics (QJ: first a transient regime is 
detected where the intensity grows exponentially, followed 
by a subsequent saturation. This regime is characterized 
by the formation of a macro-particle whose dynamics is 
responsible for the oscillations observed at the intensity 
level. Importantly, the macro-particle is shown to increase 
in size also when operating within the framework of the 
relevant self-consistent context (see Fig. EJ. In order to 
quantitatively characterize the evolution, we define a ra- 
dius and a gyroradius associated with the inner massive 
aggregate. The radius corresponds to the standard devia- 
tion of the trapped particles, namely : 

R(t) = ^((r^-rcW) 2 ),, (9) 

where rj(t) = (6j(t),pj(t)) stands for the coordinates of 
the j-th particle in phase-space, rc(t) — (rj(t))j repre- 
sents the center of mass of the macro-particle, and QjeM 
denotes the average over the subset of particles consti- 
tuting the inner core. On the other hand, the gyroradius 
characterizes the rotation of the macro-particle and reads : 

Rgyr = \/((rG(t)~m 2 )t>t aat , (10) 

where ()t>t sat denotes an averaging over time (after satu- 
ration has been reached at t sat ), and r^j = (rG(t))t>t sat - 
As it is shown in TableQlthe radius R which is almost con- 
stant in time, increases by 28% when the control is turned 
on at A = A", whereas the gyroradius is significantly re- 
duced (by nearly a factor 6). These indicators points to 
an effective stabilization of the macro-particle dynamics 
in phase-space. 

As reference to the wave, the control induces an effec- 
tive stabilization of the intensity (see Fig. I X OB . In order 
to quantify the beneficial effects of the control we mea- 
sure two quantities, namely the mean value of the in- 
tensity (I) = (I{t))t>t sat , and the average fluctuations 
Al = (\I(t) — (I)\)t>t sa f As confirmed by inspection of 
Tabled the value (I) is increased by 34% when the con- 
trol is turned on, while the oscillations have damped by 
39%. 
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Figure 9. Snapshot of the dynamics in the (#,p)-plane at 
t — 1000, when Hamiltonian (£) is simulated with N = 10000 
particles and A = —0.1321. Here the control is switched on 
at t — 300, i.e. after saturation has been achieved. The grey 
points correspond to the particles of the chaotic sea, the dark 
ones form the macro-particle [see Eq. 0]. 
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Figure 10. Intensity of the wave in the saturated regimes 
for Hamiltonian 10 for four values of the parameter : A = 0, 
-0.037, -0.0746 and -0.1321. For each value of A, the system 
starts with a negligible value for the intensity, and a monoki- 
netic beam of particles. 



4 Towards an effective stabilization of the 
intensity 

In the previous section we provided numerical evidence 
that our method enables to modify the size and position 
of the macro-particle. Here we shall take one step further 
and demonstrate that we can alter the internal structure 
of the massive agglomerate. In particular, we will show 
that the macro-particle, originally composed of invariant 
tori, can be chaotized, thus inducing an effective mixing 
while keeping the particles trapped. 

Information on the internal structure of the macro- 
particle can be gained from the residue of its central el- 
liptic point, in the test-particle picture. In particular, the 
destruction of invariant tori can be expected if this or- 
bit turns hyperbolic Thus, the residues allows one 
to trace the following two prescribed behaviours for the 
system : 

1. the chaotization of the macro-particle is controlled by 
the residue of the central periodic point (with r.n. 1), 
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Figure 11. Residues of the elliptic (upper curve) and hyper- 
bolic (lower curve) points with r.n. 7, and residue of the peri- 
odic point with r.n. 1 (bold curve) of Hamiltonian 11H With 
all the parameters set to zero except A 3 . 

2. the trapping of particles is still guaranteed by the ex- 
istence of invariant tori at the borders of the macro- 
particle (see Sec. El- 

To this aim, the residue method is implemented with a 
control term containing additional parameters. We con- 
sider the following family of Hamiltonians : 



H c 1 J9,p,t;X) = H lp (0,p,t)-2 



N 



K 

E 

k=2 



Afe cos (k9 



(ii) 

with A = (A2, A3 ... Xk), and where H lp is given by Eq. 10. 
Hence, the original system is recovered for A = 0. 

First, the macro-particle is enlarged, while keeping its 
core composed of invariant tori : the system is perturbed 
by tuning the parameter A3 around 0; the residues of the 
elliptic and hyperbolic periodic orbits with r.n. 7 are mon- 
itored until condition Q is reached at A3 « —0.0263 (see 
Fig, lilt , which in turn corresponds to restoring the in- 
variant tori. Meanwhile, the residue of the central elliptic 
periodic point (with r.n. 1) remains stable. 

The Poincare section of Hamiltonian 111 111 with A3 = A3 
and Afe^3 = 0, corroborates these predictions (not shown 
here, but similar to Fig. though the external elliptic 
islands have r.n. 6 instead of 5). 

In order to get the macro-particle chaotic, we per- 
turb the system (by making Xq different from 0) until the 
residue of the central elliptic orbit (of r.n. 1) crosses 
or 1, whereas the residues of the periodic orbits of r.n. 6 
satisfies condition©. 

Though the latter condition cannot be exactly matched 
for small values of Xq, the residue R is shown to attain its 
minimal value at A§ = —0.0201 (see Fig. ll2t . and one can 
expect a regularization of the dynamics in this region of 
phase-space. Meanwhile, the residue of the central elliptic 
orbit has crossed 1, so one might expect chaos in this do- 
main. Notice that it was not possible to track numerically 
this central periodic orbit beyond Ag. Different scenar- 
ios are compatible with this finding. Possibly the latter 
has been broken or alternatively it is associated with a 
tiny basin of attraction which prevents the multi-shooting 
Newton-Raphson to converge. 
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Figure 12. Residues of the elliptic (upper curve) and hyper- 
bolic (lower curve) points with r.n. 6, and residue of the central 
periodic point with r.n. 1 (bold curve) for Hamiltonian 11H 
with A 3 = A§ and \k^3,& — as functions of A 6 . 
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Figure 13. Poincare section of Hamiltonian 111 IP with A3 = 

-0.0263, A 6 = -0.0201 and A fe ^ 3 , 6 = 0. 

The Poincare section for Hamiltonian lilt with A3 = 
A3, A6 = Ag and Xk^3,6 = 0, depicted in Fig. Il3l displays 
a chaotic core, while the invariant tori of its borders are 
preserved, ensuring that the particles remain trapped. We 
have computed the Lyapunov exponents in the different 
regions of phase space : Inside the macro-particle it is 
equal to 0.06 whereas, in the chaotic sea it is 0.14, i.e., 
more unstable. 

In order to complete the analysis, a set of simulations 
are performed within the self-consistent framework where 
N particles interact with the wave. The additional pertur- 
bation is naturally introduced as : 

H c N ({0j, Pj }, </>,I; A) = H N ({9j, Pj },<l>,I) 

- 2 V^^ AfcCOs(fc0J + 0Xl2) 

where Hn is given by Eq. JQ. Simulations in the two 
regimes of control, (A3 = A3, A/^3 = 0) and (A3 = A§, A6 = 
Ag,Afc^ 6 = 0), confirm the prediction of the mean-field 
framework : the macro-particle has been increased and 
homogenized (see Fig.^J- 
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Figure 14. Snapshot of the dynamics in the (<9,p)-plane of 
the N = 10000 particles at t = 1000 (the control started 
at t = 300) for Hamiltonian 112t with A3 = —0.0263 and 
A6 = —0.0201. The control term is applied at t = 300, after 
saturation has occured. The grey points correspond to the par- 
ticles of the chaotic sea, the dark ones form the macro-particle. 
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-0.0263 


-0.0263 


A 6 








-0.0201 
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0.7 


1.15 
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0.17 


0.03 


0.0015 
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0.62 
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0.97 


AI 


0.73 


0.44 


0.017 



Table 2. Characteristics of the macro-particle and the wave 
for Hamiltonian 11211 with Xk^3,6 — 0. 




l , , , , , 1 

100 200 300 400 500 600 



Figure 15. Intensity of the wave in the saturated regime for 
Hamiltonian il'Il with = 0. For each value of A, the 

system starts with a negligible value for the intensity, and a 
monokinetic beam of particles. 



According to our estimates, the radius has been in- 
creased by more than 50%, while its gyroradius has been 
decreased by two orders of magnitude (see Tab.EJ. As for 
the intensity of the wave (see Fig. IT5|l , the fluctuations 
have been reduced by a factor 40, while the mean-value 
has been raised by 50%. 

In terms of topology of phase-space, it is interesting 
to observe that the regularization of the dynamics is also 
found within the framework of the original self-consistent 
picture (see Fig. llfijl , Indeed, the confinement of the trapped 
particles for Hamiltonian lfTT|l to small domains of phase- 
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Figure 16. Poincare section of Hamiltonian 11211 with A3 = 
—0.0263 and A6 = —0.0201, when the trajectory intersects the 
plane dl(t)/dt = 0. The different trajectories are represented 
by different grey levels. 
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5 Conclusion 



In this paper, we proposed a stabilization strategy in a 
wave-particle system and in particular in the a FEL set- 
ting. The technique uses an appropriate perturbation of 
the associated microscopic dynamics. We carried out our 
approach in the saturated regime of this system where 
the macro-particle, resulting from a set of invariant tori, 
plays a crucial role. In the framework of a mean-field ap- 
proximation, a linear stability analysis of selected peri- 
odic orbits (coined by a residue method) can point at the 
creation and destruction of invariant tori, allowing to re- 
shape the macro-particle through specific values of the 
parameters of the additional perturbation. Numerical sim- 
ulations performed for the original JV-body self-consistent 
picture confirm the validity of the proposed approach and 
clearly demonstrate that the method allows one to pre- 
dict a set of parameters for which a significant reduction 
of the oscillations of the intensity is found. Note that the 
residue method is quite flexible and applies to any two 
degree of freedom Hamiltonian systems. We expect that 
other families of perturbation terms would yield similar 
results. 

As a final remark, we observed that the chaotization 
of the center of the macro-particle (by a bifurcation of the 
central elliptic periodic orbit which organizes the macro- 
particle) is associated with a strong stabilization of the 
intensity. 
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